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Abstract 

According to Dirac's ideas, the vacuum consists of infinitely many vir- 
tual electrons which completely fill up the negative part of the spectrum 
of the free Dirac operator D°. In the presence of an external field, these 
virtual particles react and the vacuum becomes polarized. 

In this paper, following Chaix and Iracane (J. Phys. B, 22, 3791-3814, 
1989), we consider the Bogoliubov-Dirac-Fock model, which is derived 
from no-photon QED. The corresponding BDF-energy takes the polariza- 
tion of the vacuum into account and is bounded from below. A BDF-stable 
vacuum is defined to be a minimizer of this energy. If it exists, such a min- 
imizer is solution of a self-consistent equation. 

We show the existence of a unique minimizer of the BDF-energy in 
the presence of an external electrostatic field, by means of a fixed-point 
approach. This minimizer is interpreted as the polarized vacuum. 



1 Introduction 

The relativistic quantum theory of electrons and positrons is based on the free 
Dirac operator, which is denned by 

3 

D° = -iJ2 a k d k + (3 ■= -ict-V + (3 (1) 

k=l 

where a = (01,02,03) and 
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We follow here mainly the notation of Thaller's book [SI]. We have chosen a 
system of units such that H = c = 1, and also such that the mass m e of the 
electron is normalized to 1. 

The operator D° acts on 4-spinors, i.e. functions tp 6 7i := L 2 (M 3 ,C 4 ). It 
is self-adjoint on TC, with domain ff 1 (IR 3 ,C 4 ) and form domain # 1//2 (IR 3 , C 4 ). 
Moreover, it is defined to ensure 

(D ) 2 = -A + l. 

The spectrum of D° is (—00; — 1] U [1; 00). In what follows, the projector asso- 
ciated with the negative part of the spectrum of D° will be denoted by P°: 

P :=X(-oo;0)P°)- 

We then have 

D°P° = P°D° = -VT^ap = -p°VT^a, 
D°(l - P°) = (1 - P°)D° = vT^A(l - P°) = (1 - P°)VT^A, 

and 

n = n°_ ®n° + , 

where H°_ := P°H and H° + := (1 - P°)H. 

The fact that the spectrum of D° is not bounded from below is the source 
of many difficulties in Relativistic Quantum Mechanics. To explain why a 
free electron does not dissolve into the lower continuum, Dirac's idea |131 114j 
was to postulate that in the absence of external field, the vacuum contains 
infinitely many virtual electrons which completely fill up the negative part of the 
spectrum of D°. This Dirac Sea should be seen as an infinite Slater determinant 
S7° = ip® A • • • A if)® A • • • where (V ; ?)«>1 is an orthonormal basis of 7i°_, whose 
density matrix is precisely 

i>l 

The projector P° is often called the bare vacuum 

Let us now add an external coulomb field, created for instance by a system of 
smeared nuclei. The density of protons in this system is a nonnegative function 1 
n such that J K3 n = Z, the total number of protons in the nuclei. In our system 
of units, the external coulomb potential felt by the electrons is — a<p, where 
(p = n * tj and a is a small dimensionless coupling constant, usually called the 
Sommerfeld fine structure constant. The Dirac operator with this external field 
is 

jjotv D o _ aip _ ( 2 ) 

Dirac postulated that the charge of the bare vacuum is not measurable. In- 
deed, P° commutes with translations, so its density of charge must be constant, 

1 However, we shall not limit to nonnegative L densities n in this paper, since the model 
we want to study is able to describe the vacuum interacting with both matter and antimatter. 
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and cannot create any electric force. However, in the presence of an external 
field, the virtual electrons react, by occupying the negative energy states of a 
new Dirac operator which does not commute with translations: the vacuum is 
polarized. This polarization of the dressed vacuum, which takes the form of a 
local density of charge, is measurable in practice. 

The vacuum polarization plays a minor role in the calculation of the Lamb 
shift for the ordinary hydrogen atom (comparing to other electrodynamic phe- 
nomena), but it is important for high-Z atoms [13] an d even plays a crucial role 
in muonic atoms |18| I24j. It also explains the production of electron-positron 
pairs, observed experimentally in heavy ions collisions [Tl l4lfl 13T1 15U1 I2U] . 

When the external field is not too strong, a good approximation (called the 
Furry picture [2T|) is to define the polarized vacuum as the projector 

P a(p :=X(-oc;0)(^)- 

Note that in reality, the polarized vacuum modifies the electrostatic field, and 
the virtual electrons react to the corrected field. This remark naturally leads to 
a self-consistent equation for the dressed vacuum, and to a fixed-point iterative 
procedure for solving it. If one starts the procedure from P , the first iteration 
gives P Q< ^, and this explains why the Furry picture is a good approximation. 
But in practice, corrections to the Furry picture are necessary for high accuracy 
computations of electronic levels near heavy nuclei. These corrections can be 
interpreted as the second iteration in a Banach fixed-point algorithm (see, e.g., 
da section 8.2]). 

In physics, self-consistent equations are usually derived as Euler-Lagrange 
equations of an energy functional. It is the case, for instance, in the nonrel- 
ativistic Hartree-Fock model |4Uj . Similarly, the self-consistent equation for 
the vacuum has a variational interpretation: it is satisfied by a minimizer of 
the Bogoliubov- Dirac- Fock (BDF) energy functional. This functional was first 
introduced by Chaix and Iracane j^j, as a possible cure to the fundamental 
problems associated with standard relativistic quantum chemistry calculations. 

In these calculations, electrons near heavy nuclei are usually treated, in first 
approximation, by the Dirac- Fock model [S3], a variant of Hartree-Fock in which 
the kinetic energy operator —A/2 is replaced by the free Dirac operator D°. 
This approach gives results that are in excellent agreement with experimental 
data [HU [2^1 E3 E] • When a higher accuracy is needed, the more sophisticated 
multiconfiguration Dirac- Fock- mo del is used to take into account correlation 
effects [2S], and one can even compute the small corrections predicted by QED 
(vacuum polarization and radiative corrections), using perturbation methods. 
However, the Dirac-Fock model suffers from an important defect: the corre- 
sponding energy is not bounded from below, contrary to the nonrelativistic 
Hartree-Fock case, and this leads to important computational difficulties (see 
[Hj for a discussion and detailed references). From the mathematical viewpoint, 
one can prove that the Dirac-Fock functional has critical points which are so- 
lutions of the Dirac-Fock equations |15| 144] . but these critical points have an 
infinite Morse index, and the rigorous definition of a ground state is delicate 
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The second problem with Dirac-Fock is its physical derivation: one would 
like to interpret this model as a variational approximation of Quantum Elec- 
trodynamics (QED), which is believed to be the exact theory. An interest- 
ing attempt in this direction has been made by Mittleman |42| . but it is not 
fully convincing. According to this author, the ground state of the Dirac-Fock 
model should be obtained by means of a max-min procedure applied to the 
no-photon QED Hamiltonian Hqed- hi this procedure, the reference projec- 
tor for the normal ordering of Hqed is not fixed, and the vacuum polarization 
terms are neglected. Then one has to look for a projector which maximizes, in 
the Hartree-Fock approximation, the ground-state energy of the normal-ordered 
Hamiltonian. From a mathematical viewpoint, Mittleman's max-min principle 
has been investigated in the papers [3J I17[ Ej. In the case of zero or one 
electron [UE1> it works very well and one shows that the projector P a<f> of the 
Furry picture is the optimal reference. But it seems, from the counterexample 
given in [6], that serious problems occur when there are several electrons. 

In their work 9 , Chaix and Iracane derive their new mean-field model 
(which they call Bogoliubov-Dirac-Fock) from the no-photon QED Hamiltonian, 
normal-ordered with respect to a fixed reference: the free projector P . They 
keep the vacuum polarization terms, pointing out that they are "necessary for 
the internal consistency of the relativistic mean- field theory and should therefore 
be taken into account in proper self-consistent calculations, independently of 
the magnitude of the physical effects" 9 , page 3813]. This allows them to obtain 
a bounded-below energy: the ground states are simply defined as minimizers 
and no max-min procedure is needed. A minimizer without charge constraint, 
if it exists, is a projector satisfying a self-consistent equation: it should be 
the negative spectral projector of the mean-field Hamiltonian generated by the 
nuclear charge density, corrected by a vacuum polarization effect. This self- 
consistent projector is the stable dressed vacuum. Now, if one restricts the BDF 
functional to the charge sector — N, and if one can find a minimizer, it will be 
solution of the Dirac-Fock equations for N electrons, corrected by a vacuum 
polarization term :9,, section 4.2]. The Dirac-Fock model is thus interpreted 
by Chaix-Iracane as a nonvariational approximation of BDF. In other words, 
the Euler-Lagrange equations only differ by small terms, but the variational 
structure is completely different since the DF functional is strongly indefinite 
(i.e., it is not bounded below and all its critical points have an infinite Morse 
index). 

As we have seen, the Chaix-Iracane model has several advantages as com- 
pared to the standard Dirac-Fock model: it is more accurate (taking into ac- 
count vacuum polarization effects), its physical derivation is more convincing, 
and the ground state solutions have a simple definition as minimizers of the 
BDF functional. The drawback is that it is not easy to give a meaning to the 
quantities (energy of the vacuum, charge density of the vacuum) appearing in 
the BDF model. It is well known that there are divergent quantities in QED 
even after normal ordering, but Chaix and Iracane do not address this problem 
in their work. The first rigorous works on the BDF model are due to Chaix- 
Iracane-Lions |1U| and Bach-Barbaroux-Helffer-Siedentop ~4 a . In particular, in 
[I], the authors give a rigorous meaning to the BDF energy in the class of op- 
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erators with trace, and show that it is bounded below if one fixes the reference 
for normal ordering. Then Bach et al. vary the reference for normal ordering 
and neglect the vacuum polarization terms, which are experimentally small and 
mathematically divergent. This approximation is exactly the one made by Mit- 
tleman [32j in his formal derivation of the Dirac-Fock model. In the present 
work, our approach is different: we keep P° as reference for normal ordering, we 
study the full Chaix-Iracane model of the dressed vacuum (without neglecting 
any divergent term), and we control the divergences thanks to a momentum 
cut-off. Nevertheless, the paper jl] has been an important source of inspiration 
in our study: it contains very useful mathematical ideas and results, in partic- 
ular the lower bound on the energy (see Theorem ^ in the present paper, which 
is a mere rephrasing, in our framework, of this estimate). 

Mathematically speaking, we shall say that a vacuum is an orthogonal pro- 
jector P with the additional requirement that 

Q = P-P° e& 2 (H) , 

where 62 (W) denotes the class of Hilbert- Schmidt operators on TC. As explained 
in the Appendix, this condition guarantees that P is the (unrenormalized) den- 
sity matrix of a dressed vacuum in the electron-positron Fock space associated 
with the free projector P°. This dressed vacuum may be seen (formally) as an 
infinite Slater determinant $1 = ipx A • • • Aipi A • • ■ where (V'i)i>l is an orthonor- 
mal basis of PTC. Since the model takes the free vacuum as reference according 
to Dirac's ideas |131 114j . Q is the true (renormalized) one-body density matrix 
of f2. Following [2J (with notations from 4J, the BDF energy of the dressed 
vacuum can be written (formally) as follows, as a function of its renormalized 
density matrix : 

S(Q) = tr(D°Q) -a\p Q ^\\\ ^f^dx dy - f // ^fdx dy 

" ' U (3) 

with pq(x) = tr C 4 Q(x, x). By formal computations, Chaix and Iracane [HI 
Section 4.2] show that the Euler-Lagrange equation of this functional is 

[P,D Q ] = 0, (4) 

where 

D Q := D aLp + ap Q * -L - a^^-, (5) 
H \x-y\ 

Q = P — P°. For a minimizer, the second order condition implies a more precise 
relation between P and Dq, which takes the form of a fixed-point equation: 



P = X(-oo;0)(A?) = X(-oo;0) (d** + ap Q * ^ - «^^|) 



(6) 



Remark that if (p = (no external potential), then P° is already a solution of 
this equation since Q = in this particular case. 
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The present paper contains the first mathematical study of a fixed-point 
algorithm for finding a solution of (|Bj). Notice that the use of a fixed-point 
method to solve a self-consistent equation is very common in quantum chemistry 
and physics and that most of the numerical algorithms used in practice are based 
on this idea. For a mathematical existence result using the Schauder fixed- 
point theorem, see the resolution of the Hartree equations in [57] . See also |S] , 
where rigorous results are given on the convergence of standard Hartree-Fock 
iteration schemes. For the determination of a projector in no-photon QED, the 
fixed-point method has been used for the first time by Lieb and Siedentop |39j . 
Their goal was to replace P° by a new (self-consistent) projector commuting 
with translations, as reference for normal ordering in the absence of external 
field. We use the Banach fixed-point theorem as in pH], but our physical model 
is very different, and the necessary estimates are much more delicate in our 
case. 

Of course, we have to make an assumption on the external potential: it 
should have a certain regularity, and should not be too strong, otherwise we are 
not able to prove that the iteration method converges. If one is only interested 
in the existence of a minimizer, it is possible to remove the smallness assumption 
on the potential, but for this purpose the constructive fixed-point approach must 
be replaced by a direct - and non-constructive - minimization argument |28j . 
The regularity assumption cannot be dropped: this is a well known phenomenon 
in QED when P° is chosen as reference for normal ordering (see e.g., |35|). But 
this regularity is not really a restriction from the point of view of physics: 
point-like nuclei do not exist in nature. 

In [3], the operator D°Q is assumed to be trace class, so that the expressions 
|2J) and (JSJ) are well defined. Unfortunately, it turns out that, when ip is nonzero, 
Q = P — P° is never trace class if P is a solution of © • Therefore no minimizer 
can exist in the trace class &i(TC) in the presence of an external field. So we 
must try to define the BDF energy and the self-consistent equation for operators 
which are not trace class, and this leads to several difficulties. 

A first problem occurs with the definition of tr(D°Q) in (J3J). To solve it, we 
will have to extend the trace functional to a bigger class of compact operators, 
namely the operators with "P°-trace" (see Section \'2. II below) . 

A second problem occurs with the definition of the density pq. For this rea- 
son, we introduce a momentum cut-off A, which means we replace the ambient 
space H by 



Since D is a multiplication operator in Fourier space, Tt\ is invariant under P 
and we keep the notation P° for the restricted operator. With the cut-off, the 
integral kernel Q(x,y) of P — P° £ ©2(Wa) becomes smooth for any dressed 
vacuum P, and one can easily define Pq{x) = tr C 4 Q(x,x). Notice that, even 
with our ultraviolet cut-off, Q = P - P° is never trace class if P is a solution 
of (jHJ) and if an external potential is present. As we shall see it later on, our 
results will be valid under a technical condition of the form ay/In A < C for some 
constant C. For a small a, this leads to an extremely large A, which corresponds 
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to scales that are far beyond the reach of experimental and theoretical physics 
at the present time. But our conditions do not allow to pass to the limit of an 
infinite cut-off. 

Note that if one expands the right-hand side of equation (jSJ) in powers of the 
small parameter a, the first order term contains an expression which diverges 
logarithmically as A goes to infinity. When the exchange term Q{x,y)/\x — y\ 
is neglected, a simple algebraic manipulation allows to rewrite a posteriori our 
cut-off version of equation © in a renormalized form, with the divergent term 
removed, and the "bare" constant a in front of the charge densities replaced by 
a smaller, "dressed" , coupling constant 

a 

adr "l + |flogA 

(details will be given in a forthcoming paper 28 ). The dressed constant is the 
observable one. Its experimental value is ad r — 1/137. This kind of "charge 
renormalization" associated with a momentum cut-off is standard in the physics 
literature (see, e.g., jHHl Equation (7.18)]). With this interpretation, the limit 
case of an infinite cut-off appears as unphysical (it would correspond to Qd r = 0, 
which means no more electrostatic interaction). 

Remark 1. In the Furry picture, that is to say when P = p a <£^ it is known since 
the very beginning of QED [TIJ EHJ 1221 ESI EU that the density p aip associated 
with Q aLp = P aLp — P° is never well-defined if no ultraviolet cut-off is imposed. 
One possible regularization procedure ^jl E3 EH] is to remove the divergent 
part of p a,f , which is (formally) proportional to the nuclear charge density n. 
This gives a renormalized density p" e n which can be defined without the help of 
a high momentum cut-off. This procedure has recently been clarified by Hainzl 
and Siedentop in [22]. Some interesting features of />ren, in the case of strong 
external fields, were obtained by Hainzl in [2Zj- We do not want to give a precise 
definition of p?en here and we refer the reader to [2H1 123 • 

It would be tempting, instead of using a cut-off, to renormalize pq a priori 
in equation (jEJ), as in [2H] . But we do not know how to solve the resulting 
renormalized equation © if no cut-off is made. Moreover, even if we could find 
a solution without momentum cut-off, its interpretation as a minimizer of the 
BDF energy would be unclear. 



The paper is organized as follows. In the next section, we define the 
Bogoliubov-Dirac-Fock model and state our main results. For the sake of clarity, 
we have brought all the proofs together in Sections El and 0J In the Appendix, 
we explain in our language, for the reader's convenience, how the BDF energy 
is deduced from no-photon QED by Chaix-Iracane in 

Acknowledgment: C.H. wishes to thank Heinz Siedentop for suggesting him 
the possibility of studying a self- consistent model of the polarized vacuum, dur- 
ing his time as post-doc at the LMU (Munich). The authors acknowledge sup- 
port from the European Union's IHP network Analysis and Quantum HPRN- 
CT-2002-00277. E.S. acknowledges support from the Institut Universitaire de 
France. 
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2 Model and main results 



In this section, we study the Bogoliubov-Dirac-Fock model introduced in [Qjlllj. 
Our system of notation is similar to 0] , with the difference that we keep all the 
terms describing the vacuum polarization. This forces us to deal with operators 
which are not trace class, unlike [I]. 

2.1 An extension of the trace functional 

In order to give a meaning to the expression "tr(D°Q)" even when Q is not 
trace-class, we need the notion of "P°-trace". In this section only, we work in 
an abstract Hilbert space I). 

Definition 1. Let P be a projector such that P and 1 — P have infinite rank, 
and A £ 62(f))- We shall say that A is P-trace class if and only if A ++ := 

(1 — P)A(1 — P) and A := PAP are trace class. Then we define the P-trace 

of A by 

tip(A) := tr(A++) +tr(A__). 

We denote by <Sf(f)) the space of all Hilbert- Schmidt operators which are P- 
trace class. 

Notice that if A is a trace class operator, then A £ ©f (h) and to(A) = 
trp(A) for any projector P. 

Remark 2. In |54[ Section 5.7.2], a similar definition in connection with su- 
persymmetry is made and the name "supertrace" is used. 

The following result, whose proof is given in Sectional will be used repeat- 
edly in the sequel. 

Lemma 1. Let P and P' be two projectors such that P — P' € ©2(f)). Then 
A is P-trace class if and only if it is P' '-trace class, and in this case top (A) = 
top, (A). 

Another useful fact is that when A is Hilbert-Schmidt and A + P is a pro- 
jector, then A has a P-trace, as explained below: 

Lemma 2. Let P and P' be two projectors on a Hilbert space, such that P'—P is 
a Hilbert-Schmidt operator. Then P'—P is P-trace class. Moreover, top{P' — P) 
is an integer which satisfies 

to P {P' - P) = to ((P' - P) 2n+l ) 

for all n > I, and to P (P' - P) = when \\P' - P|| 6oo < 1. 

In our framework, a consequence is that, for any vacuum P such that Q = 
P — P° G 62 (Wa), trpo(Q) is an integer which can be interpreted as the charge 
of the dressed vacuum P (see the Appendix for comments in this direction). 
When P solves the self-consistent equation © and ip is not too strong, we will 
see that P is close to P°, so that its charge will be zero, according to the lemma. 
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2.2 The Bogoliubov-Dirac-Fock model 

As in [3], we are going to extend the BDF energy to a convex set of compact 
operators, which can be interpreted as one-particle density matrices of quasi- 
free states. This kind of extension is standard for mean-field models depending 
only on the one-body density matrix (see [3H1 El E] ) ■ 

In the whole paper, we assume that the nuclear charge density n = —Atp/An 
belongs to the Hilbert space 

£={/ GL 2 (IR 3 ,1R), £>(/,/)< 00} , 

where 

f(k)g(k) 



£>(/,<?) =4tt J 



C| 2 

We will choose the following Hilbert norm on <£ 



-dk. 



1 + W 2 \f(k)\ 2 dk^ 2 



\k\ 2 

The Bogoliubov-Dirac-Fock energy is defined by 

£(T) = trpo(D°r) - aD( Pr ,n) + -Dfapr) - % [[ ^'^ dxdy (7) 

2 2JJ \x-y\ 

on the set 

Qa ■■= {r e &(°(H A ) I - p° < r < 1 - p°, Pr e <t} . (8) 

in O, 

Mk) = ^ TV C4 (f (p + k/2,p - k/2)) dp 

is the Fourier transform of the charge density pr, which, formally, is the diagonal 
of r € 62 (Wa), as explained in the Introduction. Thanks to the momentum 
cut-off, pr is compactly supported, so that pr £ L , hence 

pr(x) = Tr c4 (T(x,x)) = — [ f Tr C 4 (f(p,q))e ix ^dpdq. (9) 



(2tt) 3 



HM<a 



Clearly, the function T £ 62 (H-a) l— *• Pr 6 Cq n L 2 (R 3 ) is continuous. Notice 
that if for instance f,g£ i/ 1 (]R 3 , M), then the electrostatic energy is simply 



D(f,g)= If 



^^ dxdy, 
\x - y\ 



but for functions in (£, it does not necessarily have a meaning as a Lebesgue 
double integral in direct space. 

Note that the set Qa is convex, and that the elements of P° + Qa are not 
necessarily projectors. In fact, it is an easy exercise to show that an element of 
Qa is extremal if and only if it is of the form P — P°, with P a projector. It 
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then follows from Lemma [2 that the set of all extremal points of Ga coincides 
with Q a =Va- P°, where 

V A = {P orth. projector | Q = P - P° G ©2 (Ha), Pq G £} • 

It will turn out that, under some assumptions on a, A and n, the BDF functional 
has a unique minimizer on Ga which is extremal. As a consequence, 

mf{£(P - p°), Per A } = inf{£(r), r g Ga} ■ 

In the next subsection, we give necessary and sufficient conditions satisfied 
by a minimizer of £. 

2.3 Study of the BDF energy 

We first state the following result, which is an easy translation, in our frame- 
work, of the stability estimate proved by Bach et al [3] (see also jlUj): 

Theorem 1. Let be n G £. Then 

1. £ is well-defined on Q A ; 

2. ifO<a<-, then 

Oi 

VrG£? A , £{T) + -D(n,n)>0 (10) 

and therefore £ is bounded from, below on G\, independently of A; 

3. ifO < a < ~ and n = 0, then £ is non-negative on Ga ^} \1(JI/. being the 
unique minimizer. 

Remark 3. Note that the result is optimal in the sense that the functional 
becomes unbounded from below when n = if a > 4/-7T, as shown in ^U] and 

Remark 4. Since aD(n,n)/2 in (J10[) is the electrostatic energy of the field 
created by n, ()1U|) means that the total energy of the system is nonnegative. 

Proof of Theorem^ We only explain here why £ is well defined on Ga, the rest 
being identical to the proof of Theorem 1 in [3] (see |U Eq. (18)-(19)]). 
If r G SC'^'Ha), then we have 

P°D°TP° = D°P°rP° = D°T ++ G &i(Ha) 

since P° commutes with D° and \D°\ < Vl + A 2 . With a similar argument for 
1 - P°, we obtain that D°T G Sf°(W A ). Therefore, tr P o(£>°r) is well-defined 
and 

tr P o(Z) r) = tr(D°r ++ ) + tr(£>°r__) = tr(|,D |r++) " tr(|Z>°|r__) (11) 
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(notice that, due to the constraint — P° < F < 1 — P , one has T ++ > and 
r < 0). On the other hand we have by Kato's inequality 

ff lT{X ^dxdy<^tr(\D^) 
J J \x-y\ 2 

showing that this last term is well-defined since \D °| is bounded on TCa and 
r e 6 2 (H A ). □ 

We are interested in minimizers of the BDF functional, and we expect them 
to be in the class Q\ = Va — P°- This leads to the following definition 

Definition 2. We say that a projector P is a BDF-stable vacuum if and only 
if P — P° is a minimizer of £ on Q\. 

When there is no external potential, P° is the unique BDF-stable vacuum 
|10| Hj. which corresponds to Dirac's ideas. But if we consider a non- vanishing 
external potential (p = n * A, then P° obviously cannot be BDF-stable, since 
it is easy to create a state — P° < 7 < 1 — P° such that £(7) < = £(0). This 
means that the vacuum is necessarily polarized. 

More precisely, one can easily derive necessary conditions satisfied by a 
BDF-stable vacuum P. To this end, a perturbation of the form Q + 7 = 
P — P° + 7, with 7 E &i(7i\) such that — P < 7 < 1 — P is considered in 
Chaix-Iracane formula (4.8)], and the energy £(Q + 7) is expanded to get 

£(Q + 7)=tr(D Q 7) + |D(p 7 ,p 7 )-| JJ h ^ y) J dxdy + £(Q), (12) 

a formula which is valid when 7 G Si (Ha)) the operator Dq being defined in 
©• 

Remark 5. In jlj Formula (21)] and [7], the polarization potentials appearing 
in Dq and the energy of the vacuum £(Q) were neglected by the authors who 
used the following functional 

£p( 7 ) = tr (D«* 7 ) + |£>(p 7 , P7) " f // T^rf^ rf f > ( 13 ) 

with the constraints 7 G Si and —P < 7 < 1 — P (and even ^7(1 — P) = 
in [7]). Then a procedure taking the form supp inf_p< 7 <i_p £p(7), related to 
Mittleman's work |32], was considered in ,4 a . For the case of the vacuum (no 
constraint on the trace of 7), the solution is the Furry picture P = P atp with 
7 = 0, as shown in [3J. We refer the reader to O page 3809] and ^7J Ej 
for comments and results concerning Mittleman's max-min in the case of N 
electrons (which corresponds to the additional constraint tr(7) = N). 

^,From formula ()12|) . it can be seen that a BDF-stable vacuum must satisfy 
the fixed-point equation ©. The converse is also true under some assumptions: 
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Theorem 2 (BDF-Stability). Let be P G Va and n E €. We assume that 
there exists a positive constant d such that 

d\D Q \ > \D°\ with ad^ < 1, (14) 

where Dq is defined in @ . Then, the following assertions are equivalent 

1. P fulfills the equation 

P = X(-oo;0)(A?) = X(-oo;0) +a P Q*y- { - « ■ (15) 

2. P is the unique BDF-stable vacuum, i.e. P - P° is the unique global 
minimizer of £ on Qa- 

The proof of this theorem will be given in Section 3. Some arguments are 
directly inspired of [I]. 

Remark 6. One can try to go further: if we translate the ideas of Chaix and 
Iracane in our language, the ground state of a molecule consisting of nuclei with 
total charge density n, surrounded by a cloud of N electrons, should solve the 
following constrained minimization problem 

min{£:(r), r e Ga, tr p0 (r) = N}, 

with JV£N\{0}. If this minimization problem has a solution Q, it will solve 
a self-consistent equation of the form 

Q = X(-^)(d q )-p°, 

where \i G (— 1; 1) is a Lagrange multiplier associated with the charge constraint, 
and interpreted as a chemical potential. For a not too strong external field 
(p, it should be possible to prove that \x £ (0; 1) and that the vacuum II = 
X(-oo;0)(-^q) stays neutral, which means trpo(II — P ) = 0. Therefore, we 
could split P = Q + P° = X(-oo;/i)(Dq) i n the form 

N 

p = n + £|Vfc)Mfe|. 
fc=i 

The mono-electronic wave functions ipk would be solutions of the Dirac-Fock 
equations (with high momentum cut-off), perturbed by vacuum polarization 
terms: 

D Q ip k = e k ip k , < e k < 1 , l<A:<iV. 

It is our goal to study this constrained variational problem in the near future. 
The present work, which deals with the unconstrained case, is a first step in 
this direction. 
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2.4 Existence of a BDF-stable vacuum 



We may now state our main Theorem. We recall that the norm on (£ is 



2 \ 1/2 

-\f(k)\ 2 dk 



i±!2!_im.^ 

I*l s 

Theorem 3 (Existence of a BDF-stable vacuum). Zei be n € £ and 

6 G (0; 1). TTien for all A and a such that 

2^y^ra\\n\\£ < b and a < at(A) (16) 

where 



a b (A) 



C(l - b) 



°° v^gA' 

t/iere exists a unique BDF-stable vacuum P, which is a solution of 

( 1 Q(x y)\ 

P = X(-oo;0) (P>Q) = X(-oo;0) ( -° + «(PQ ~ n ) * ~ | ~ " ^ _'^| J (17) 

mf/i Q = P — P° . Moreover, we have trpo(Q) = 0. 

Remark 7. The first constraint 2y / 7ra||n|| l r < b means that the external field is 
not too strong. It explains why a neutral polarized vacuum is obtained (since 
trpo(Q) =0). In our proof, this constraint on the external field is necessary for 
the fixed point algorithm to converge. The second constraint a < ab(A), which 
essentially reduces to a^/\ogX < C(l — b), is a technical condition due to our 
choice of norms, but we were unable to drop it. It disappears if the exchange 
term ^ J J ^i|z^r dx dy is neglected in the energy, as can be seen from the proof. 
A precise definition of the constant C appearing in this result is given in the 
proof. 

Remark 8. There is an interesting symmetry property of the solutions of (|17jl 
when n is replaced by —n. Namely, if P is a solution of (|17f) with external 
density n, then P' = Q' + P° is a solution of (|17|) with external density —n, 
where Q' = —CQC' 1 , C being the charge conjugation operator 54, page 14]. 
The two dressed vacua P and P' have the same BDF energies and satisfy 
Pqi = ~pq, as suggested by the intuition. For this symmetry between matter 
and antimatter to be true, it is essential to have the Fermi level at and not 
at —1 (see, e.g., the comments of jlHl page 197] about this fact). 

2.5 Idea of the proof of Theorem OH the fixed-point algorithm 

We end this section with a brief description of our fixed-point algorithm, used in 
the proof of Theorem 13 and which could be useful for practical computations. 

A natural scheme for solving (fTTj) would be to construct a sequence (Qj)j>o C 
Qa by taking Qq = and 



Qj+i = X(-oo;0) (d° + a( PQj - n) * 1- - «^?^) " P°- 



(18) 
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Expanding this expression in powers of a, and considering the total density 
one can write the following recursion formula in Fourier space: 

oo 

P^Q~+i(k) = —n(k) — aBx{k)p'q. {k) + a p~ifi(Qj)(k) + a n PrJ&hPq?) {k) ■ (19) 

n=2 

The notations B\ , p\$ and p n are defined precisely in Section 14.11 (see the 
subsections 14.1.31 and l4*7l.4j) . The important point is that B\(k) is a positive 
function which diverges logarithmically as A —* oo for any fixed k, whereas the 
other terms stay bounded. From (fT9|) we thus see that the scheme (fT8|) would 
converge under a condition of the form a log A < C. 

To improve this condition, we use a better algorithm in our proof. Our 
modified scheme consists in defining a sequence of pairs (Qj, p'j)j>o c5a x ^ 
such that (Qo,p'o) = (0, —n) and 

< ' Q ^ = *<-«*» { D ° + a ^ * n " a i^j ) " p0 (20) 

, p'j+i = L 0> a )Pq j+1 + C 1 - L {ot,A))p'j 

where p'q.(x) := p Q \x) — n(x) = tr C 4 Qj(x,x) — n(x), and L(q,A) is the linear 
operator which, in the Fourier domain, is just the multiplication by the function 
(l+al?A(&)) -1 • The second equation in the iteration scheme (|2fl|) can be written 
in the form 

oo 

p^(fc) = (1 + aB A (k))- l [-n(k) + a p^Q 3 ){k) + ^a" Prfh^j>(k)] (21) 

n=2 

The divergent term now only appears in the denominator. So one expects a 
much better convergence. In the proof of Theorem we show the convergence 
of the algorithm (|2()|) under the conditions (jlfij) but we believe that it converges 
independently of the cut-off A. 

It can be seen from our proof that this holds when the exchange term 
a -x^~r i s neglected in (|2Uj) . In this case, the algorithm converges independently 
of A to the solution of a reduced fixed-point problem (without exchange term) , 
which is the unique minimizer on Q\ of the convex functional 

£ red (r) = trpo(D r) - aD(p r ,n) + ^D{p r ,p r ) . 

3 Proof of Theorem [2] 

In this section, we prove Theorem|21 To this end, we first need to prove Lemmas 
□ and El 
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3.1 Proof of Lemmas [T] and [21 

Proof of Lemma^\ Let be P and P' two projectors such that P' — P £ ©2(f)), 
and a Hilbert-Schmidt operator A which is P-trace class. This means that 
PAP and (1 - P)A(1 - P) are trace class. 

Let us first show that P'AP' is trace class. To this end, we write 

P'AP' = (P' - P + P)A(P' - P + P) 

= (P' - P)A(P' - P) + (P' - P)AP + PA(P' - P) + PAP. (22) 

This shows that P'AP' is trace class since the last term is in ©1 by assumption, 
P' — P and A are in 62 j and P is bounded. The same computation shows that 
(1 - P')A(1 - P') is trace class. 
We now compute 

tr[P'AP'] = tr[(P' - P)A{P' - P)\ + tr[(P' - P)AP\ + trfPy^P' - P)] + tr[PAP] 
= \x[A ((P' - P)(P' - P) + P(P' - P) + (P' - P)P)] + tr[PAP] 
= trL4(P'-P)]+tr[PAP], 

where we have used the formula tr(AB) = tr(BA), valid for A,B E 62- The 
same computation gives 

tr[P' + AP' + ] = tr[A(P| - P+)] + tr[P + AP + ] = -tr[A(P' - P)} + ti[P + AP+], 

where we have used the notation P + = 1 — P and P' + = 1 — P' . Summing this 
two results, we obtain the formula trpL4] = trp/L4]. □ 

Proof of Lemma\B We introduce B = P'-P. We have B 2 = P' — P'P - PP' + 
P = (1 - P)P(1 - P) - PBP. This implies that (1 - P)P(1 - P) and -PBP 
are non-negative trace class operators. We now use the proof of 2, Theorem 
4.1]. Since B £ ©2, we infer P 3 £ ©1 and so (P', P) is a Fredholm pair, in the 
language of [2]. Therefore, tr(P 3 ) is an integer and satisfies tr(P 3 ) = tr(P 2n+1 ) 
for all n > 1. Now we have 

P 3 = P 2 P' - P 2 P = P'BP' + PPP 

Applying this result to 1 — P' and 1 — P, we find 

P 3 = (1 - P')P(1 - p>) + (1 - P)P(1 - P). 

Summing this two identities, we obtain by Lemma ^ 

2tr(P 3 ) = tr P /(P) +tr P (P) = 2tr P (P). 

This shows that trp(P) indeed equals the index of the pair of projectors (P', P) 
defined in [2j, an integer which vanishes when ||P' — P|| e < 1 by the results 
of 0. □ 
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3.2 Preliminaries 

To prove Theorem I2J we also need the following 

Lemma 3. Assume that tp = p * ry for some p E (£. T/ien 

\\ H1 = Mp\\ € , \<p\\ LOB < C^Me, M L e < C 6 47T\\p\\ € 



where : = 27r i/2 anc ^ ^6 is the Sobolev constant for the inequality \\u\\ L 6^3^ < 

^6 ||VM|| L 2( K 3). 

Proof. We have 



l + |fc| 2 .~, ,_> 



i,,|2 I dfc = tot / 1*1 (! + 1*1)1 wi 



and so 



(27r) 3/2 ll^ll^- (2vr) 3/2 VJ R 3 |^(1+|^ 



1*» - 7^372- I^Hl 1 ^ 7^372- L mrrrn^ ) 47r IWI< 



The rest is easily obtained by the Sobolev inequalities. □ 

Lemma 4. Let P be a projector in V\ and Q = P — P° . Then Dq is bounded. 

Proof. Due to the cut-off in Fourier space, D° is bounded on 7i\. On the other 
hand, if p = p* rr for some p E £, then 99 £ L°° by LemmaEland so this is also 

a bounded operator. Let us now denote R(x,y) = ^-y\ • ^ e then have 



Q(x,y)f{y) 



dy 



\x-y 

and since, by Kato's inequality, 

2 



IQ(^y)l 2 ^.7 / \f(y)\ 2 



F - S/l / Wr3 F - 2/| 



/ l»%<f</,l^l/> 

Jr3 \x-y\ 2 



lQ{X,y) } 2 dxdy< -tr(\D°\Q 2 ) (23) 
\x — y\ 2 

this shows that R< C \ D°\ 1 / 2 and so R is bounded. □ 

Lemma 5. Let P be a projector in V A and Q = P-P°. Then D Q T E ©f°(HA) 
for all r E £/a <fid we /lave 

tr P o(£)°r) + a£>(pg - rc,p r ) - atr (^3^^ = tr P o(r> Q r). 

Proof. Remark that i?r = is trace class, since , L and , r ^',F L are 

in S 2 by (El . Let us now define L> = D~ a{pQ ~ n) *V\ = D° + a(p Q - n) * X 
and P' = X(-oo ; o)(-D)- By the result of Klaus-Scharf 35] (see also [HOI and the 
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proof of Theorem OJ), it is known that P' - P° G 6 2 (W A ). Thus P'DTP' = 
DP'TP' G ©i(Wa) since r G 6f = ©f by Lemma □ and L> is bounded by 
the proof of Lemma |IJ Therefore, DqT = DT + RT is in &^°. 
To show the expected equality, we prove 

tr P o(L>°r) + aD(p' Q ,p r ) = ti pQ (DT) (24) 

where p'q = pg — n G C This will end the proof since the other term is trace 
class. The general idea of the proof is to approximate T by a trace class operator 
for which this equality is true, and to pass to the limit. However, the behaviour 
of the associated density in the space £ is not obvious and to overcome this 
difficulty, we shall also approximate the density p'q to obtain a potential in 
L 2 (R 3 ). We thus start by choosing a sequence pj which converges as j — > +00 
to p'q in £, such that cpj = pj * rr is in L 2 (R 3 ). We can choose for instance 

Pj( k ) = PQ^^m^i/j)- We now show 

tTpo(D°T)+a p r ipj =tr P o(D j r) (25) 

for all r G Q\, and where Dj = D° — acpj. To this end, we may find a sequence 
r™_ of finite rank operator which converges to = (1 — P°)TP° in ©2- Then 

r n •= ( r++ r +-^\ 

• \{Ti_y r__; 

converges to T in ©2- Since r n G Si for all n > 0, we have 

tr p0 (L>°r n ) + a / pr"^' = trpo^T"). (26) 

./R 3 

By Q, the function Q G 62(Wa) ^ PQ ^ L 2 (M 3 ) is continuous. Therefore, 
pr™ - * Pr in L 2 (R 3 ). Since tpj G L 2 (IR 3 ), we may now pass to the limit in (|26(1 
and obtain 

lim ( tr P o(L>°r n ) + a / p r ™¥>i ) = tr P o (D°T) + a / pr</?j, 
n ^°° V Jr 3 / ./r 3 

where we have used that 

trpo(L>°r n ) = tr(P°r™ + )+tr(L>°r™_) = tr(D°r ++ )+tr(D°r__) = tr P o(L>°r). 

Let us now pass to the limit in the right hand side. Indeed, we can write, by 
Lemma ^ 

trpoODjT") = tr^CDj-r*) = tr^Pjl^Pj) + tr(Dj(l - Pj)r n (l - Pj)) 

where Pj = X(-oo;0)(-^i) an d since Pj -P° G 6 2 by [Mi- Now, using (|22|). it is 
easily seen that Pjr n Pj -» PjrPj and (1 - Pj)r n (l - Pj) -» (1 - Pj)r(l - Pj) 
in Si as n — > 00, since this terms can be expanded as a sum of trace class 
operators and products of at least two Hilbert-Schmidt operators converging 
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strongly in ©2- Since Dj is bounded by the proof of Lemma 0J we obtain that 
ti P i(DjT n ) -^n^oo tTpi(DjT) = trpo(-DjT) by Lemma^ 

As a conclusion, we have proved ()25j) for all F £ Q\. To finish the proof, it 
remains to pass to the limit as j — ► +00. Since 




prifj = D(p r ,pj) 



and pr £ £ (recall that T € £/a.)j Pi — * p'q strongly in £ as j — > 00, we may pass 
to the limit in the left hand side of (|25[). To pass to the limit in the right hand 
side, we use again the fact that 

trpopjT) = tTp'XDjT) = trpjPjTPj) + tv(Dj{l - Pj)T(l - Pj)). 

By the results of Klaus-Scharf jSEj (see also the proof of TheoremEJ, it is known 
that Pj — P' — > in ©2, since pj — ► p'q in £. Using again (|22|). it is then easily 
seen that PjTPj -> PTP' and (1 - P-)T(1 - Pj) -»■ (1 - P')T(1 - P') in Si as 
j — ► 00. Since Pa,- — > P in ©oo by Lemma @J we may thus pass to the limit and 
obtain the desired equality (|24jl . □ 

3.3 End of the proof of Theorem [2] 

We start by proving 1) => 2). We thus consider a projector P that satisfies 
the assumption of the Theorem, and is also a solution to the equation P = 
X(-oo;0)(Mg)- We fix some F £ Q\ and show that £ (F) > £{Q). To this end, 
we write £ (F) = £(Q + F') where F' = F — Q = F + P° — P. By assumption, 
T fulfills -P° < r < 1 - P°, and so F' fulfills -P < T < 1 - P. Using Lemma 
El we may expand £(Q + F') and obtain 

£(Q + r') = trpo (D Q F') + ^D{p T ,,p r ,) -^jf lT ^ y J dxdy + £(Q). 

Using now LemmaQ] we see that it is thus sufficient to prove that 

trp (D Q F') + ^D(p r ,p r ) -^JJ ^yj dxdy > 0, 

for any F' £ ©f (Ha) such that F' / 0, p r > £ £ and -P < F' < 1 - P, which is 
an easy adaptation of the proof of [IJ Theorem 2]. 

We now show 2) 1). Let be P which satisfies the assumption of the 
Theorem, and such that Q = P — P° is a minimizer of £ in B\. We therefore 
have, by formula (|12jl . 

tr(P Q7 ) + ^D(p y ,p 7 ) -^ff ^^J dxdy > (27) 

for all 7 £ ©1 (TC\) such that — P < 7 < 1 — P. The proof of Theorem 4 by Bach 
et al. 4 now implies that P = X(-oo;0){Dq)- Their proof is done with D aip 
instead of Dq but they also mention that it can be extended to a more general 
case, provided ^ ct(Pq) and P, 1 — P leave the domain of Dq invariant, which 
is the case here. □ 
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4 Proof of Theorem [3] 



In this section, we prove Theorem |3] by using a Banach fixed-point method. 
4.1 Preliminaries 

We start by defining the norms and spaces that will be used to apply this well 
known result. In fact, one of the main difficulties we faced in this work consisted 
in finding suitable Banach spaces. 

4.1.1 Norms and spaces 

We choose the following norms 

Q-=(JJ E(p-Q) 2 E(p + q)\Q(p,q)\ 2 dpdq ) j ' , 
\R\n -=(JJ E J: P {p ~ € \R(P, g)| 2 dpdq ' 



+ 



2 \ 1/2 

-\p(k)\ 2 dk 



E ( k ) ,cm2 



where 



|j :r : | / \k\*E(k)*\m\ 2 dk' "" 



E(k) = ^/T+W 



and denote by Q, 1Z, <£ and y the associated Hilbert spaces. The dual space £' 
of £ will be also useful and we introduce 

2 \ 1/2 

\((k)\ 2 dk 



E(k) 

In the following, it will be easier to use the norm H-Rli? where R = ®( x ' y } 
°> " 11 /v \x—y\ 

in our estimates and a relation with ||Q||q will then be needed. To this end, 
we first need the following well known Lemma, which will be useful throughout 
the rest of the proof. 

Lemma 6. For all £ and r] in M 3 , we have 

Vs > 0, E(Z) S < 2 5(s) (£(£ - 7]) s + E{ v ) s ) (28) 



Vs G R, E(£) s < 2^E(C - r,) s E{ri) 



\\»\ 



(29) 



with 



5(s) 



s if < s < 1 

s - 1 if s > 1 
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Remark 9. A trivial consequence of (|28(1 is the following inequality 

' <nhn(— — r,-=7-^ r J . (30) 



E(p) + E(q)- \E(p + q) } E(p-q) 



We shall also need the following 
Lemma 7. We have 



E(p + g) 
,gGR3 E(p) 2 E(p - qf 



Proof. Let us introduce the function f(p,q) = E(p + q)E(p) E(p — q) for 
(p,g) £l 3 xK 3 . We have 

o < /(p. ») = = g " < 2 * (2 *l E(p - $ < 1 



E{p) 2 E{p - q) 2 ~ E(p) 2 E(p - q) 2 ~ E(p)E(p - q) 

by Lemma El Therefore lim(p i9 )_KX3 f{p, q) = and / attains its maximum 
on M 3 x R 3 . Computing V g /(p, q), we see that at a critical point of /, p 
and q are always parallel. It therefore suffices to study the function g(x, y) = 
E(x + y)E(x)~ 2 E(x — y)~ 2 for (x,y) E RxR. It is then easy to see that 
max R 2 g < 2 (the bound (|3Tj) is indeed not optimal). □ 

Now we can give a connection between \\R\\n and ||Q||q when R = 
(we also recall the easy relation between \p\<t and \\ip\\y when ip = p* A). 

Lemma 8. If p G £ and Q £ Q, then we have ip = p * yjy G 3^ and R(x,y) = 
^ ^ anc ^ more precisely 

\\<f\\y = M\p\k, 

\\R\\n < C R \\Q\\ Q , (32) 
||^| J Dor 1 || e2 < V2||i?|k, (33) 

with 

C R := -X; inf sup f £(2x) e / — — — ^ = ) . (34) 

Proof. We have 

ft, x i / Q(p-i,g-i) 
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so we obtain, for some fixed 9 £ (0; 2) 



112 ff E (P 



\2 



= 8 JJ ^^-\R(u + v,u-v)\ 2 dudv 

8 /■/■/•/• £(2u) 2 Q(Z + u, Z - u) • Q(Z' + u, Z' - u) 



(2tt 2 ) 2 7777 E(2«) \l-u\ 2 \l'-u\ 



8 /■/■ £(2?;) 2 
' ' dudv — ; — - x 



du dv dl dl' 



(2tt 2 ) 2 7 7 #(2u) 

E{2l)^Q(l + v,l-v) E{2l')^Q{l' + v,l'- v) , 
E(2l')^\l - u\ \V - u\ E(2l) 1 ^ L \l - u\ \V - u\ 

< » //// w if c + »■; - f 

~ (2vr 2 ) 2 7777 #(2u) E(2Z') 1+e |Z-u| 2 |Z '-«| 2 

< 8 // £(2u) 2 £(2Z)|Q(Z + «, Z - w)| 2 ^(0 dt> dZ 



where 

/'y (2Z)" /" /" 1 

^(0 -=-^2 J J E (2u)E(2l>y+e\l - u\*\l> - u\* dudl ' ■ 
Now, let us introduce 

dit 



C e := sup ( £(2x) e / 



£(2u) 1+e |u- jg| 



Remark that 



du 1 f du 

< 



E{2u) l + e \u- x\ 2 ~ 2 1 + s \x\ e 7 R3 H^lu-exl 2 
where e x := showing that Cq < oo when 6 (0; 2). Now we have 

E(2V)® f 1 / f 1 

= p / dn ^(2n)^|Z-.| 2 ^(2Z0^|Z'-n| 2 

E(2Z) e /• 1 /^V 

" Wf J E(2uV+e\l - u\* X Ce ~ \2^J 



dl' 



and so 

\\Rfn < 8 (^) 2 / /" i?(2^) 2 i?(2Z)|Q(Z + v, I - v)\ 2 dv dl < Q^Q " \\Q\\ 2 Q 

which ends the proof of (|32|) . 

To prove (|33l) . we remark that we have, by (|31l) . 

II pin l-ill 2 ff \R(p,q)\\ , ff E{p-q) 2 \R( P ,q)W 
W R \ D o\ \\e 2 =JJ^P^ d P d ^ 2 JJ E(J+q) dpdq - 
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4.1.2 An estimate from below for Dq 

We now state a Lemma in which we give a lower estimate for the operator 



:= D + afj, * -j— r 



Dn a := D + an * r-r - a 

\r\ \x — y\ 

by D°, in terms of the spaces introduced above. 

For our result, we are interested in Dq = DQ tPQ ^ n but this definition with 
an arbitrary density fi will be useful later on. 

Lemma 9. Assume that (Q,[a) £ Q x <£ are such that 

a(2^Mc + V2C R \\Q\\ Q ) < 1. (35) 
Then Dq^ is a bounded operator which satisfies 

\D Q J > (l - a (2v^H|c + V2C R \\Q\\ Q )) \D°\. (36) 
Proof. We have, with (// = /j,* A, 

l^'wl^ < Iv'ILoo HIl 2 ^ 271-1/2 HI e II I A) I -u\\ L 2 

by Lemma |21 and 

\\Ru\\ L 2 = {{RlDol-^DoluWv < || J R|Z)or 1 || e2 HAiMI^ < V2C R \Q\\ Q \\\D \ ■ u\ L2 

by (E3J). This shows that \<p' - R\ < (27r 1 / 2 || / u|| £ + V2C r \\Q\\q)\D°\, the square 
root being monotone. This proves that Dq^ is bounded since Dq is bounded 
on Tt\, and gives the expected inequality. □ 

Remark that Lemma |§] will be useful when we shall apply Theorem [2] (see 
the condition ()14|) in the statement). It also implies ^ o~{Dq^), a fact that 
will be used to compute the projection X(-oo ; o)(A2,/J- 

4.1.3 Expansion by Cauchy's formula 

We want to solve the equation 

Q = X(-oo;0) (d° + a(p Q - n) * -L - 4^j) - X(-oo;0)(D°) := F 1 (Q). 

If a (2y/n\pQ - n\<i + y/2C R \Q\o) < 1, then $ ct(Dq) by LemmaEl We 
may thus use the method of and expand F\ by Cauchy's formula 

-i r+CG / -i 1 \ °° 

27r 7-oo \ D Q+*V D u +irjJ ^ 



n=l 



where 



1 f + °° If 1 \" 
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pQ = pQ- n , 1 pq = pq*m' RQ{x,y) = ■— — ; 



We shall write 



Qr 



A;,// k+l=n 



JUJ={l,...,n}, 
|/|=fc, |J|=i 



+00 



+ ir? y V D° + iriJ 



where Rj = Rq if j £ I and i?j = (p'q if j £ J (Qk,i is the sum of all the terms 
containing k Rq's and I <£q's). We also denote := pQ kr 
Hence our equation can be written 



Q = J2 an Qn(Q,P Q ) 

n=l 

00 

PQ = J2a n pn(Q,P Q ), 



(37) 



n=l 



where we recall that Q n and p n depend on both Q and p'q = pQ — n. In order to 
have a better condition on a and A, we shall now change the second equation 
for the density, by taking into account the special form of the first order term 
p\. To this end, we need to compute this term explicitely. 



4.1.4 The first order density 

Recall that 



Qo,l — 7T 



+00 



dr] 



1 



1 



2vr 1 D° + ir? y D° + irj 



so that 



r+00 ± ^ 

Q ,l(p,q) = (2tt)- 5 / 2 / dr? ——<p Q (p-q) — — . 

7_oo a-p + (5 + irj ^ a - q + + ir) 



We now introduce 

M(p,q) 



dr] 

00 

^ r+00 



1 



1 



a ■ p + f3 + it] a ■ q + (3 + it] 
a • p + (5 — irj a • q + (3 — irj 



-00 
1 



p 2 + 1 + rf 



+ 1 + r] 2 



(a-p + P){a-q + (3) ^ 



Hence 



E{p) + E{q) V £(p) 



QoAp^q) = 2 5/ 27r 3/2 ^qCp- q) m (p^q)- 



(38) 



23 



This enables us to compute 

m(k) = £^Tt c «(QiTi(i + fcAz-fc/2))di 

1 



|1[<A 

J£(ife) / Tr C 4 (M(l + k/2,l- k/2))dl 



2 



16vr 3 

-7 Q {k)B A (k) (39) 



where 



d m L_ /" (/ + k/2) • (Z - fc/2) + 1 - g(j + fc/2)£(Z - fc/2) 

A vr 2 |fc| 2 i| Z |< A £(Z + fc/2)£(Z - k/2)(E(l + fc/2) + E(l - k/2)) ' 

(40) 

This function is computed in jl^j 

1 z 2 - z 4 /3 



B A (fc) = - / 

Jo 



l-z 2 l + |£f(l-z 2 )/4 
and it is logarithmically divergent since 

Ba(0) = - 4^ dz = T~ lo s( A ) " #■ + T" lo § 2 + (VA 2 ). 

7T Jo 1 — Z z 37T 97T 3"7T 

4.1.5 Equation 

We are now able to introduce the function on which we shall apply the fixed- 
point theorem. According to what we said above, the equation in pq can be 
written 

oo 

P Q = Ys an Pn(Q,p' Q ) (41) 
n=l 

or equivalently (we forget the dependence in Q and p'g for simplicity) 

oo 

PQ(k) = -aB\(k)p' Q (k) + ap{fl(k) + a n p^(k) 

n=2 

and 

— 1 1 / °° \ 

= - l + aB A (kf {k) + TT^aW +E«W) J ■ ( 42 ) 

which is more adapted to a fixed-point argument, since the divergent term 
B\(k) appears now in the denominator. 

Notice that we now study an equation for the full density p'q = pq — n and 
not for pq as previously. We therefore introduce the following space: 

X = Q x € 
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consisting of all the pairs (Q, p') such that Q 6 Q and // G £. Notice that in 
this space, p' can be different from pq — n. However, we shall find a solution of 
the equations in this space, which satisfies p' = pq — n. We also introduce on 
X the norm 

\\(Q,p')\\=C R V2\\Q\\ Q + 2^\\pk, 

where we recall that Cr is defined in Lemma 03 In the following, we shall keep 
the notation p' to remind the reader that the equation indeed concerns p'q and 
not pq. 

We now introduce the function F : X — > X defined by 
F(Q,p') = (fq(Q, P ') , F p (Q,p')j 

where 

oo 

Fq(Q, P ') = X(-oo ; o)(Aj,pO " ^° = J>"Q n (Q,p') (43) 

n=l 

. I if 00 \ 

KM = -TT^ S + TT^T, + E«"S«3,p') . (44) 

\ n=2 / 

Qn(Q,p') and p n (Q,p') being defined in (|3*7|) (replace p'q by p'). Remark that 
Pn = Pq„ for all n > 2. In the proof of Theorem |31 we solve the fixed-point 
equation in X 

F(Q,p') = (Q,p'). 

4.2 Existence of a fixed-point of F 

To prove our main Theorem, we need the following estimates 

Proposition 10. Assume that (Q,p ! ) £ X is such that ^ o~(Dq p /). Then we 
have 

+00 

\\F(Q, p')\\ < 20F||n||ct + «i(A)a \\(Q, p')\\ + \\(Q, p')\\) n (45) 

n=2 



\F'(Q,p')\\ < Ki(A)a + a^TiK„ (a || (Q, p')||) n 1 (46) 



aK n (a ||^,p / )|| vn 

n=2 

where 



/ a \ f^W? /-— /- , VlogA\ CW2 /-— 

«i(A) = max I — ^=— vlog A , \/2C R + ^-j= I ~ A ^oo ^_ V lo g A 

and (K n )n>2 a sequence of positive numbers independent of A and which 
satisfies K n ~ n _ +ao Ky/n for some constant K. 



25 



To prove this proposition, we have to do some tedious estimates. Before 
starting this proof, let us show that Theorem |3] follows from Proposition I1UI 

Proof of Theorem^ We introduce the function f(x) = Y^=2 K nX n , which is a 
power series with a radius of convergence equal to 1. The estimates 1)45(1 and 
(|46|) can be written 

\\F(Q,p')\\ < 2^F||n|| £ + Kl (A)a \\(Q,p')\\ + f {a \\(Q,p')\\) (47) 

\\F'(Q,p')\\ < Ki(A)a + a/' (a \\(Q,p')\\). 

To apply the Banach fixed-point theorem, we now have to find a ball B(0,R) C 
X which is invariant under the function F and on which F is a contraction. Let 
R > be some fixed radius. We have 

sup \\F'(Q,p')\\ < Ki(A)a + af'(aR) := p. 

(Q,p')eB(0,R) 

Moreover, we also have 

\\F(Q,p')\\ < \\F(Q,p')-F(0,0)\\ + ||F(0,0)|| < p \\(Q,p')\\ + \\F(0,0)\\ ■ 

Therefore a condition for the ball B(0,R) to be invariant under the action of 
F is |F(0,0)|| < (1 - fJt)R. Notice that since F(0,0) + (0,0), this inequality 
also contains the contraction condition p < 1. Additionally due to Lemma 
we assume aR < 1 as well as 

an 



4(1 - aR) 

due to Theorem |21 and Lemma El which is equivalent to 

1 

a < ; . 

- tt/4 + R 
As a conclusion, if (a,R) fulfills 

2V7r|n|<£ + olRki(A) + aRf'(aR) < R 
1 

7r/4 + i?' 



1 (48) 
a < 



then we are able to apply the Banach Fixed- Point Theorem on B(0, R). Remark 
that these inequalities also contain the conditions p < 1 and aR < 1. Notice 
also that if (a,R) is a solution to (|48jl. then (a',R) is a solution to (jjHj) for all 
a' < a, since the function which appears on the left of (|48|) is increasing in a. 
Now, if we assume that 2y / 7ro!|n|£ < b, we obtain that if (a, R) fulfills 

- + aRn\ (A) + aRf'{aR) < R 

a i (49) 

a ~ vr/4 + i?' 
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then it also fulfills (|48|) . The first inequation of (|49[) is simpler when it is written 
in terms of the variables a and x := aR. It becomes 



- + K 1 (A)x + xf'(x) < - (50) 
a a 

which implies x G [b; 1). Now, given 6, A and x let us call ab t \(x) the maximal 
value of a such that (|5U|) holds, i.e. 

ab < A(x) = ^(Afr+t/'fr) ' (51) 

which is defined for x in [b;l). Since \im x ^i ^^(x) = a^A^) = 0, we may 
denote by x\ £ (b; 1) the largest maximizer of the function ab a m the interval 

We now define Rb(A) := x\/ ab^{x\) and 



o*(A) ==min(^( 3 : A ), ?r/4+ 1 iZ6(A) ) 



As a conclusion, for all < a < a;,(A), (a, i?;,(A)) is a solution of (@HJ). This 
means that F is a contraction on B(0, Rb(A)), on which we can apply Banach 
Theorem. This gives a unique solution to the equation F(Q,p') = (Q,p') in 
B(0,R b (A)) CX. 

Let us now show that P is indeed a solution to ©. In fact p' is a solution 
to (|42|) and so p = p' + n is a solution to (|41j) . On the other hand, we have 
Q = X(-oo;0) (Agy) — -P ! and (|4"T|) means exactly that p = pq. Hence, P is a 
solution to P = X(-oo;0)(Dq)- Thanks to the proof, we know that P satisfies 
the assumptions of Theorem |2 and so P is the unique BDF-stable vacuum (i.e. 
P — P° is the unique global minimizer of the BDF energy) . 

To end the proof, let us study the behaviour of a&(A) as A — > oo. Computing 
da^A^/dx, we find that x\ must satisfy the equation 

K 1 {h) + fM= Xk{X l~ h) f"(x h ). 

Since ki(A) diverges as A — ► oo, we see that f"(x\) — > oo and therefore xa — > 
1 as A — > oo. Now, since /'(x) = (f"(x)), we obtain that /'(xa) = 

o A ->oo(«l(A)). Thus 

QVH^a) ~A^oo — ttt and Rb{A) ~a- 



ki(A) UK ' " 00 1 



As a conclusion, 



1-6 v^F(l - 6) C(l - b) 

" fe(A) «i(A) C^v^ " "^gX 

withC = ^. □ 
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4.3 Proof of Proposition [TOl estimates 

In this section, we prove the claimed estimates of Proposition 1101 We will have 
to introduce many constants. For the sake of clarity, a guide is provided to the 
reader at the very end of the proof, section 14.3,51 
Remark first that we have 

1 < 
1 + aB A (k) ~ 

for all k G R 3 . Therefore, to estimate the norm \\F p (Q, p')\\(r, it suffices to 
estimate the norms of pifl(Q,p') and p n (Q,p'), due to (fl4"|) . 

For (Q, p') £ X, we introduce the notation R = £ 72, and ip' = p' * A- 

We then remark that 

\\F(Q,p')\\ <2v^M c + al(Qa,pi )0 )| + 5> n \\(Qn,Pn)\\ 

n>2 

and estimate each term separately. A similar argument can be done for F'(Q, p'). 



4.3.1 First order terms 

Lemma 11. We have the following estimates: 

\\Qo,i\\Q < {l0g ^ l,2 y\y = 2(logA) 1 /2|^| k! 

|Qi,ols < \Rq\h < c R \\Q\\ Q , \\ Plfi k < Cr{1o ^ A)1/2 \\Q\\q- 

Therefore 

||(Qi,Pi,o)|<«i(A)||(Qy)||, 

where 

«i(A) = max (^0oiA, ^2C R + . (52) 

Proof. Recall that 

Qo,l(P> a ) = 2 5/2 3/2 ^'(P ~ l) M (P> ( 53 ) 

where the matrix M(p,q) is defined in (|38[). and whose properties are summa- 
rized in the following 

Lemma 12. Let A+(p) = and A"(p) = ~ (< " P 2 +^ £(p) 6e ifce pro- 

jections matrices in C 4 onto £/ie eigenspaces of D° in Fourier space. We then 
have 

TV C4 (M(p,q)) = ~± E{p) l + E{q) ^ (A + (p)A-(?)) 
\M(p,q)\ 2 = Tr c4 (M(p, q)M(p, q)*) = 8 i Tr C 4 (A+(p)A-(g)) 
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TV C 4 (A+(p)A-(g)) = Tr C 4 (A-(p)A+(g)) =1 ' ' 7 + ' 



£7(p)S(g) 
Moreover, we have 

Vp, g G M 3 , Tr c4 (A+(p)A-(g)) < min ( ^ ( fc~^ /2)a , 2) . (54) 

Proof of Lemma WA We only prove (|54l) . We have Tr C 4 (A + (p)A~ (g)) < 
|A + (p)| |A~(g)| = 2, so when i := gwfe^jT^p — 1> there is nothing to prove. 
Now we have 

1 _ pj_g_+l 1 _ I 2 - fc 2 + 1 ^ = y + g fe = p-g 



E{p)E{q) E(l + k)E{l - k) 

1 



V(l + t) 2 -4zt 



where t = \k\ 2 /E(l) 2 and z = n$n^ G [0; 1). When i G [0; 1) and 2 G [0; 1), 
the expression above is decreasing in z and so we obtain 

p-q + 1 l-t 2t 

1 - — < 1 , = < 2t 

£(p)£(g) " v/flTO 5 1 + 1 ~ 

which ends the proof. □ 
• Let us now treat Qo,i- From (|53j). we obtain 

|Q^(p,g)| 2 = ^3|^| 2 (p-g)|M(p,g)| 2 , 

and so 

II Eip-qfEip + q^^q^dpdq 

dkduE(k) 2 E(2u)\^'\ 2 {k) x (\u\ < A)\M(u + k/2,u- k/2)\ 2 



1 

2V 



<tJ^( I dk\k\ 2 E{k) 2 W\ 2 (k)^ • 1 ' 



2 5 vr 3 V7 W ' V 7 \ v y H <A J E(2 U ) J E( U ) 2 y / 

by Lemma IT2*1 Now we have 

/ = 4^ ( ^argsh(2A) + -Urgth ( ^ ) ) < 2vr log A 

J\u\<A E{2u)E(uy \2 V3 \a/1 + 4\/A/ / 

(55) 

for A > 3. So we obtain 

IQo,i|Q<^^|^| y = 2(logA) 1 /2| p '| £ . 

• pifi and Qifi. We have 

1 /• +0 ° 1 1 

Qi,0 = -7T- / rf 1 nn I ■ 



2tt J.oo 'D° + i7] D° + ir] 
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so that 

Qi,o(p,q) 



1 



/+oo 
dr, 
-oo 

1 1 

~2E{p)+E{q) V E(p) 



-R(p, q)- 



1 



a ■ p + /3 + ir] a. ■ q + (3 + ir/ 



and 



\Ql,o(P,Q)\ 2 < rE^ ^7^ l%>g)| 2 < -F7r^2l%^)| 2 , 



showing that 
Now, we have 



1 



(2vr) 3 / 2 
1 



25/2^3/2 



so we obtain 



WQxAq < \\R\\n < C R \\Q\\ Q . 



TT C 4(Qi !0 (l + ^,l-~}yi 

Try (R(l + ^l-^\ M (l + ^,l-^) A),//. 



I^WI< ^372 (j[ 3 ^(20 _1 |^ + fc/2,/-fc/2)| 2 d/) 7 x 



E(2l)\M(l + k/2,l - k/2)\ 2 dl 



1/2 



and finally 



^o(*)| a dfc<5B-8l^ 

Z J 7T d 



1 



<A £?(20^(/) 2 2%^™ 



log Vn 2 



by (|55|) which implies 



|pi,ok< t- \Q\q- □ 



4.3.2 Second order terms 

To simplify the presentation, we introduce the following notation: 



S p>g := (4vr)(27r) 



-3/p 



E(u)i 



1 



+oo 



(56) 



Let us recall the following inequality |48| Theorem 4.1] 

ll/(^)5(-iV)|| ep < (2vr)- 3 /f |/| LP(R3) 

which implies 



(57) 



< 



p.ap 

4tt 



II /II 



LP 
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On the other hand, we shall often use the following trick 

(E(p) 2 +rj 2 )(E(q) 2 +rj 2 ) = (E(p)E(q)) 2 + (E(p) 2 + E{q) 2 )r, 2 + r, 

> \E(p + q) 2 + ^E(p + q) 2 V 2 

> ^E(p + q) 2 E( V ) 2 



by Lemma El This implies 



1 2 

< 



v^Cp) 2 + vW E (i) 2 + v 2 ~ e (p + i) E ^) ' 



Recall now that we have Q2 = Q2,o + Qi 1 + Qo,2 with 
1 f +OD 111 

Q2,0 = -7T- / ^ n n , .._ -R „n , 



2it J_ OQ D° + irj D° + i7] D° + in 
1 Z" 1 " 00 111 



2vr J_ 00 'D°+ir] D° + it] D° + irj 

1 /'+ 00 ,1,1 1 



2vr 7_ 00 ' D° + in^ D° + in D° + £77 



2tt./-oo 1 D° + i^ D° + i^ D° + in' 
We shall now treat each term separately. 
Lemma 13. VFe /mue po,2 = and the following estimates: 

\\Q2,o\\q < 2 5 / 2 K 3/2 (C R ) 2 \\Q\\ 2 Q , iQulg < 45 6 C 6 ir3 /2 CH||Q||Q||p'| 



||Qo, 2 ||q < 2v / 10S 6 ,5CmC 6 ||p /|12 



., ., S6^6(^i?) 11^1,2 11 11 ^ .S&^CmC^Cr.. ,.. .. .. 

\\p2fl c < K? q> Pi,i e < 4 — \\P KW 2i 

7T 7T 



and so 



||(Q2,p 2 )|<«2 ||(Q,pO| 2 , 



with 

^2 = Cq 2 CrV2 + 2^C P2 , 



Cq > = max ^ ^ 3/2 ' ' ^ J 



C p2 — ITlcLX 



where Cm is a constant defined in Lemma\14 
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Proof. Step 1 : Estimates on the exchange term Q 2 . 

• Q2,o- To estimate Q2,o, we write 



1 



and so by 



IQ2,o(p,g)| < 



23/2 /-+oo 



drj 



\R(p,Pl)\ \R(pi,q)\ 



2n E(rj)WE(p + q)V 2 J R s E(p + Pl )V 2 E( Pl + q)^ dPl 



which implies 
E(p-q)E(p + q) 1 / 2 \Q^ (p,q)\ < 2 5 / 2 K 3/2 



E(p-pi)\R(p,Pi)\ 
Eip + pi) 1 / 2 

E( Pl -q)\R( Pu q)\ 
E(pi + q) 1 / 2 



dpi 



and finally 



||Q 2 ,o|| s < 2 5 / 2 K 3/2 ||i?|| < 2 5 / 2 K 3/2 (C R ) 2 \\Q\\ 2 Q . 



• Qn- We treat for instance 

1 



+oc 1 1 1 

drj nn — —R-p^r ; — "V 3 ' 

-oo 



(2tt) 5 / 2 y.oo D° + irj^D° + {7]^ D° + ir] 
and use the same method to obtain 



£(p-9)£(p + g) 1/2 IQi>^)l ^ 
£(p - pi) 



+oo 



This means that 



|fl(p,Pi)| 



dq 

(2tt)^ 2 J. O0 ^) 

S(P! - g) 



x 



(E{ P1 ) 2 + r?yl*(E{q) 2 + ^2)1/4 

00 dr/ 



- g)| dpi. 



r! 



i 



r/ 



1 



(|D | 2 + ?? 2 ) 1/4 (|A>| 2 + T7 2 ) 1/4 



6 2 



where we have introduced R' and / defined by 

E(p- 



But now 



R'- 



R'(p,q) :-- 



E{p + q) 1 / 2 



\R(p,q)\, f(k):=E(k)\<p>(k)\. 



(iDoP+r? 2 )! (\D \ 2 + V 2 )* 



< \\R'\ 



&2 
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< \R\ 



1Z 



I / \ 

(\D \ 2 + v 2 )^ (|A)| 2 + ?? 2 )^ 
1 T / 1 



s 



(|A,| 2 + ?? 2 ^ (|D | 2 + r? 2 )i 
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If we now use inequality (|57|) , we obtain 
1 . 1 



(|I) p +r? 2)J (|£> |2 + ,/2)T 



< 



©6 



(| J D |2 + r? 2 )i 
< (2-)- 1/2 ( / T 



rl/l V2 



©12 

du 



+ \u\ 2 + rj 2 ) 3 



1/6 



5 6 



Finally since 

l/L6<c 6 ||v/|| L2 = c 6 |^|y, 

and = (47r)||/3 / || £ , we obtain 



(59) 



IIQi.ils 



< 



4S 6 C 6 
2ir 



+oo 



dry 



-ll^ll 



and 



, E{r,) 3 / 2 

|Ql J l|Q<8^6^3/2Ci ? |g| Q |p / | (I 



mp lie 



• Qq,2- Unfortunately, the method used above cannot be applied to Qo,2- in 
this case, we have to calculate this term explicitely. We can write 

Qo,2 = ^2 Ql). 

ei,£2,e3S{±} 

where for instance (by a residuum formula) 



1,2 



Q+^>,g) = (27T)- 3 



Qo^2 + (P' ^) = <3 ,2 (p> 9) = 0, 

A+ (p) (p - pa) A- (pi ) (pi - q)A~ (q) 



dpi 
dpi 



(E(p) + E(q))(E(p) + E( Pl )) 

A+(p)y? / (p -pi)A"(pi)y / (pi - g)A+(g) 
(£?(p) + E(pi))(E(g)+E(pi)) 



and similar formulas for the other Qq 1 ^ 3 ■ We now treat for instance Qq~ 2 • 
Using (|3Uj). we may obtain 

£(p-pi) 

dpi— TTTT^X 



£(p - q)E(p + q) 1 ' 2 1 Q+J- 1 < 2(2ir)~ 

x A + (p)99'(p-pi)A _ (pi 



E(p + P1) 2 / 3 

£(pi -g^foi -g)| 

S(pi) 1 /3 jE ( g )l/2 



dpi. 



So, we may write 

HQ{riis<2 

where 



|L>o| 1 / 3 '' \D \V 2 



< 



©2 



2lMy| 



6, 



1 r* 1 



|D |V3 J |D |l/2 



6, 



(60) 
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Lemma 14. When Mt is defined by formula (|6'fl|) . then 



11 /"e 2 - 47r 



\k\ 2 \m\ 2 dk 



1/2 



^ere C M := 2 (/ °° fi(2 ff fi(t)i 



1/2 



~ 2.1589. 

Proof. We have |A+(p)A~ (q)\ z = Tr C 4 (A + (p)A~(g)) and so, by (J5IJ), 

|p-g| 2 |J(p-g)| 2 



dpdg|M/(p,g)|^ < (2tt) 



-3 



dpdq 



2E((p + q)/2) 2 E(p + q) 4 /z 



< (2vr)- d / dA^I/WI' 



< (2tt)- 



Finally, since by (|59*)) 
1 „ 1 



A)| 1/!j/ |A)| 1/2 



< 



6> 



dk\k\ 2 \f(k)\- 



H/3 J I n„U/ 



2£(2w) 4 / 3 £(w) 2 
t 2 dt 



E(2ty/ZE{t) 2 ' 



□ 



|A)| 1/3 |A)| 1/2 



< 



©6 



^6,5 
47T 



IIL6 > 



we obtain 



II Qo^2 Is - 2 'S'6,5CmC 6 ||p / |||. 
This result is immediately extended to the others terms and since we can prove 

\Q\ 2 = \Q++- + Q+—f + |Q-++ + Q + | 2 + |Q+-+p + |Q-+-| 2 , 
we arrive at 

/ii 2 



||Qo,2||q < 2v / 105 6 ,5CmC 6 || / 9 / ||^ 

Step 2 : Estimates on the density p2- Let us now treat the density p2- 
The general idea of the proof is to estimate (p2, m terms of the norm ||C||c' 
by using 



Kp,C)| = |Tr(QC)l= TrlQC 



Jr3 v / 



< 



\Q((p,p)\dp. 



This can be done if we know that QQ G <&\. But we have 



lock 



Q\Dq 



6i 



< \\Q\ D o\ 2 \\ 62 



IAiP C 



6 2 

<£;(A) 2 ||Q|| e2 ^||CL 2 , 



showing that Q( £ Si when ( € I 2 . So, in what follows, we shall assume that 
C 6 £' n L 2 and prove a bound depending only on || C II £' - By the density of 
<£' n L 2 in this will give us a bound on \\p\\^. 
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Let us remark first that po,2 vanishes. Indeed we have 
1 



and 



(2tt) 3 / 2 



Ip[<A 



Tr C 4 [Qo, 2 {p + k/2, P -k/2) )dp 



T^C 4 (Qo,2(p,q) 



1 
1 

x — - — 



cfy / dpiTrj 



1 



A)(p) + ^ 



p'(p-pi)x 



ip'(pi - q) 



(27T) 



+oo 



c?r/ / dpi 



x TV 



c 4 



(A)0) - ir)))(D (pi) - irj)(D (q) - irj) 



Now the terms linear in the Dirac matrices are traceless and the remaining 
terms are odd in r\ and vanish after integration. This can be easily generalized 
to po 2k f° r all fc, and is known as Furry 's Theorem in the physics literature 

123- ' 

• P2,o- We use here a method similar to what we have done above. We estimate 
for some (GC'nL 2 and Q c := Q 2 ,oC 



\Q c (p,p)\ < (2tt) 



-5/2 



dr] 



\R(p,Pi)\ \R(pi,P2)\ \((P2-P)\ dpi dp 2 
yjEipf + vWHPi) 2 + VWE(P2) 2 + V 2 ' 



< 4(2?r) 



-5/2 



+00 



dr/ 
Eft) 



dpi dp; 



E(p - pi)\R{p,pi)\ Ejpi - P2)\R(pi,P2)\ 
E{p + p l ) l l 2 E{p l +p 2 ) 1 / 2 



\C(P2 ~P)\ 



E(p 2 - p)(E(p) 2 + V 2 )^(E(p 2 ) 2 + ^2)1/4 ' 



< 4(2?r) 



-5/2 



+oo 



drj 
Etf) 



dpi dp. 



R'(p 1 Pi)R'(pi,p 2 )C'(p 2 -p) 
(E(p) 2 + rf) l / A (E{p 2 ) 2 + rf) 1/A 



where R'(p,q) = g %7jSff }l and ('(k) = E(ky 1 ((k). This means that 



\(P2,0,C)\ < 



4 f + °° dr] 
2^ 



-oc E(V) 



R'R' 



(|Do| 2 + »7 2 ) 1/4 (|D |2 + ^2)1/4 



Si 



< 



45, 



6 



(27r)(4vr) 
by l|59|). showing that 



+oo 



dr] 
E(r]) 3 / 2 



\ R '\\e 2 \W\\l* - 



SqC%{Cr) 2 K z /2 



7T 



IQIIICk 



,, I, S 6 C e (C R ) 2 2 

||P2,0k S 3 A 3 /2|y||Q. 



7T 
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• pi i- Unfortunately, as for Qq )2 , we have to calculate pn explicitely. Let us 
start for instance with pf±~, the density associated with one of the two terms 
of Qi i 



(2^ 3 / 2 / 



dpi 



A + (p)R(p-pi)A (piV(pi-g)A (g) 



(£(p) + £(g))(£(p)+£(pi)) 
We use the same method as above and estimate for some ( € <£' n L 2 the term 

A+ (p)fl(p, pi)A~ (pi)^(pi - p 2 )A- (P2) 



Q f (p,p) = (2vr)- 3 
by 

IQc(p,p)I < 



dpi dp 2 



(£7(p)+£7(p2))(£;(p)+£(pi)) 



"C(P2 ~P), 



\A+(p)R(p, Pl )A-(pi)\ W{pi-p 2 )\ 
E(p + Pl y/ 2 E{p l Y/ 2 E{p 2 ) 1 /^ 

\Q(jp 2 -p)\ x \A-{p 2 )A+{p)\ 



E(p + p 2 ) 2 / 3 



-dpi dp 2 , 



< 



0(0 x-a/2 ([ E{p-pi)\A+(p)R{p,pi)A-{pi) 

2(2vr) // — s(p+pi)^ — 

with /(A;) := £?(fc)|^(fc)| and C'(jfe) := |C(fc)|/-E(*0- Now 
\Q ( (p,p)\ <2(2^)- 3 / 2 



#1 (P> Pi ) -^TZ^TT^ \ i /a M C (P2 1 P) d Pl d P2 



with 

Ri(p,pi) ■-- 

We thus have 



E{piy/ 2 E{p 2 y/* 

E(p-pi)\A + (p)R{p,pi)A-{ Pl 
Eip + pi) 1 / 2 



< 2||A + J RA-|| 7 e||M ( ,|| 62 

< 2||A + J RA-|| K ||M c || e2 
Sq^CmCq 



6 1 

■J: ; 



|A)I 1/2 |A)| 1/3 

11/2 ^ I 



< 



2tt 



A 1 1/2 |A)| 1/3 
\\A+RA-\\ n \\C\\e\\p'k, 



6, 



6 6 



and finally 



Sq$CmC§ IIA+ if ,., 
2^ |A + i?A l^lplk- 



We now treat and estimate 

A+(p)R(p,pi)A-(pi)£i(pi -p 2 )A+(p 2 ) 



Q c (p,p) = (2tt)- 3 JJ dpidp 2 



(E(p) + E(pi))(E(pi) + E(p 2 )) 
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C(P2 -P), 



by 

/(Pi -P2)|A~(P1)A+(P2)| C'fe ~p) 

£(pi + P2) 2 / 3 ^fe) 5 / 6 
Using the same argument as above, we arrive at 

11 +-+11 ^ Sq^CmCq .. , _ .. .. ,., 
Iff ^ — ^ \\A + RA \\ n \\p 1^. 

To treat pfi~ , we remark that 

1 1 

< 



dpi dp 2 



(E(p) + E(p 2 ))(E( Pl ) + £(p 2 )) " (£(p) + E( Pl ))(E( Pl ) + £(p 2 )) 

1 

+ 



(E( P ) + E( Pl ))(£;(p)+E(p 2 )) 

and use the same estimates as above to get 

11 ++-11 s ^CmCg + _ , 
\\Pii h< A^-RA p k- 

' 7T 

Finally, since Z)ei,e 2 6{±} |A ei -RA £2 |^ = \\R\\k, we end up with 

\\ PlA \k<2^^\\p'\\J I! W^R^ln) < ^ CMC6CR \\ P 'k\\Q\\ Q - 

IT \ z — ' / 7T 

\61,62e{±} / 

4.3.3 The general n th order case 

Now that we have explained how the proof works for the second order, let us 
estimate the general n th order term. 

Lemma 15. We have the following estimates 

Vn>3, \\Qn\\Q<nKzC Q \\(Q,p')\\ n , with C Q = V2 (^^J \ 

Vn>5, \Pnh<nKn±iC p \\(Q,p')\\ n , with C p = (f^) ' 

IMk < C P4 ||(Q,P')| 4 » with C Pi '■= K2S ^ C& {1^=^ ■ 

Therefore, 

Vn>4, \\(Qn,Pn)\\<Kn\\(Q,p')\\ n 

with 

k 4 = AC R K 2 C Q V2 + 2^C P4 , K n = nC R KnC Q V2 + 2nKn+iC p ^. 
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Remark that it can be proved that K, 
behaviour for n n as n —* oo. 
Proof. Step 1 : Estimates on the exchange term Q n . 
• Qk,i with k > 1 and + / = n > 3. Recall that 



-n^oo which gives the claimed 



Qk,, 



1) 



i+1 



2tt 



JUJ={l,...,n}, |J|=fc, \J\-- 



+oo 



drj 



D° 



where i?j = Rii j £ I and i?j = <// if j G J. For the sake of simplicity, we treat 
only 



Qk,i- 2^f 



+oo i 

dr\- 

oo 



R- 



1 



D° + ir] \ D° + ii] 



1 



D° + 77? 



We have 

|8£>,9)I< 
i 



(27T) 1 - 



+oo 



31 

" 2 J — OO 



r/r/ 



fc-i 

n 



(f;(p)2+ 7,2)1/2 

:l-R(Pi,Pi+i)l 



I-R(p,pi)|x 



(E(p i+1 )2 + 772)1/4 



n-2 

»n 

j=k 

X 



(^(p,) 2 + r ? 2 ) 1 /4 



(£(p„_l) 2 +7,2)1/4 



\(f'(p n -i ~ q)\ 



(E( Pj+1 )* + 7,2)1/4 
1 



(£(g) 2 +77 2 ) 1 / 2 



dpi • • ■ dp n _] 



so by l[5B|). 

£(p-g)£(p+g) 1/2 |^(p,g)| < 



11 j_«...M/2 11 



fc+i 

2~ £(p 

(27T) 1+ 



dr, 




-|y'(Pn-i - g)|- 



W(pj-Pj+i)\ 



1 



(Efe +1 ) 2 + 7, 2 )V4 
dpi ••• dpn-l- 



Now if we use the easy generalization of Q28JI 5 

E(p - g) < £7(p - pi) + E(pi - p 2 ) + • • • + £(p„_ 2 - Pn-l) + #(Pn-l 
we obtain by a similar argument as before 



II Q\ 



kl\\ 



< 



k+l 

2— 
~2V 



fc(C /i ) fc (47rC 00 )' 



E{rj) l+ — 



+i(c R ) k {^c 00 ) l - l s & cAj c 



dr) 



E(v) 



lQ\ k Q ¥le 
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To obtain this result, we have estimated each term containing a <p' by using 
1 , 1 



(\D \i + E( V y)V^ (| J D | 2 + £(r?) 2 ) 1 / 4 



and when E(pj — Pj+i) appears in front of a <p'(pj — Pj+i) (i-e. when j > k), 
by using 



1 



1 



(IDoP + W) 1 /^ (|Z) |2 + ^)2)l/4 



< 



7: 



(|A,| 2 + £fa) 2 ) 1/4 (| A) |2 + ^)2)l/4 



< 



6 () 



II L6 



So we have 



k+l 

< 2~ 

k + l 



IIQIIIIIP'II^ 



< 2 — n(C R ) k (4irC O0 ) l K l+ k max ( 1 
which implies 



IIQmIIq < (2) 2—71(^)^(4^)'^, max 1 



iiqi&iiaV 



Qo,n Wi/i n > 3. Recall that 



Qo,n = — ^ / dq- 



2tt 



so that 



S(p-g)S(p+ 9 ) 1 /2|Q _(p j9 )| < 



,1+3? 



+oo 



rir/ 



2 J-OO 



(2tt; 

Ip'(p-Pl)l 




dpi ■ --dpn-xX 



n-2 



n 



(£(p) 2 + r/ 2 ) 1 /4( J E( ;pi )2 + ^2)1/4 
1 1 



\ \(E{Pj) 2 + r/ 2 ) 1 / 4 



(S(p n _ 1 ) 2 + r/ 2 )!/ 4 

We now use (|59j) to bound for some /i, /2 and /3 



W(Pn-l ~ q)\ 



(E(q) 2 +T] 2 ) 1 / 4 ' 



n((|A| 2 +r/2)l/4 / - 7 '(| J D |2 +r/ 2 ) l/4 



< 



& 2 



S(7/) 3 /2(4^)3 LI 



J Hi 6 



(61) 
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to obtain 



s/2 ( r 



or 



|Qo,n|c < nK n _ 1 ^(5 6 C 6 ) 3 (47rC 00 )"- 3 ||p , ||^. 
Finally, we can write for instance (recall that Coo = l/(2y/n)) 

\\Q n \\ Q < nK n/2 C Q (C R V2\\Q\\ Q + 2^11^11^)" (62) 

where 

^ [k ( 1 SqCq ( SqCq \ 3 \ /^{SeCe) 3 
Cn = V2max 1, , = V2 n ' , 

v ^ 47TCOO V 47rC ooy / 8vr 3 /2 

and since K n j 2 > if n -i when n > 2. 
Step 2 : Estimates on the density p ra . 

• Pfc,z wii/i k > 2 and n > 3. As before we treat for instance the density p' fc t 
of the ; where the k i?'s are on the left and the I <p n s are on the right. For 
some fixed (G C'n L 2 , we introduce := Q' k We thus estimate 

1 / 1 ^ 1 \ 

X (E( Pi r+ ^)l/4 ^ ( (£ (pj) 2 + ^2)1/4 l^(Pi'Pi+l)l ( S (p. +1 )2 + ,,2)1/4 J X 

x S ((WF 11 "'^ " Pj+l)l (^FTW^J x 

We now use as before 

\C(Pn ~P)\< lCi f n ~ p) l (E(p -p 1 ) + E( P1 - p 2 ) + • • • + £(p n _! - p n )) 

E{Pn ~ P) 



to obtain 



Ifc4,i.0l <n^^(C R ) fe (4vrC oo y (i- jT ^ ) IQfcl^lCfc 



and so 



bMlk<«©^(^) fc (47rC 0O )^ + M i IQ|| fc e |p'| z £ . 
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Pi : i with I > 2. We may treat for instance with the same notation as before 
\QcM\ ^ —TUT f dr, [■■■[ — 1 \Rfa Pl )\ x 



(2^) 1 +¥7_oo 'J J (S(p) 2 +r?2)l/4 

1 n_1 / 1 1 \ 

X (i?(p 1 ) 2 +r ? 2)i/4 If ((^^jVl^'fe "^'^(^T^J X 

We now use (|61[1 and obtain 

1^,01 <n 2l ^ C,6)3 Ch(47rC eo )'- a ^ J^ -^—^ j \\Q\\ Q \\p'\\U\Ck' 



and so 

l m ,\\*<n( n A ^—^ 
• po i with I > 5. We want to estimate 



X (£( Pl )2 +?? 2)l/4 Q ((£7(^)2 + 7,2)1/4 "^' +l)l (^(p. +1 )2 + ^1/4) 

Since there are at least 6 functions, we may use (|61j) twice and obtain 



IW„0l<.»^(4^)«(i/ 



OC 



E(rj) 



— 2 I iip'IkllClk' 



and so 



\\ P 0,l\k < n^^(47rC 00 )'- 5 ^_ 2 |p'||^. 



47T 

Now since K n —i < K( n +i)/2 when n > 5, we obtain 



with 



llPnlk < nK^C p (C R V2\Q\\ Q + 2vWl<t) (63) 

For p4, we notice that po,4 = for the same reason as po,2, and that > 
Therefore we obtain 

|p 4 |k < C P4 (cW2|Q|| e + 2yfr\p'\\X (64) 
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with 



a 



Pi 



4K 2 S 6 C 6 
4ir 



max 



/ CqSq V 



7T 



□ 



4.3.4 The third order density p 3 
Lemma 16. We have 



INk<C P3 (CflV2||Q|| s + 2V^||p / |k 



and therefore 
with 



\\(Q3,p 3 )\\ <« 3 ||(Q,pO 



«3 



3CrK 3 / 2 CqV2 + 2^C P3 , C P3 



15C M 56(56,4) 2 (g 6 ) 4 



Proof. Notice that thanks to the previous proof, we already have some estimates 
on p3 t o, p2,i and pi ;2 - It remains to study /9o,3- As before and as in |2H], we 
have to compute po^ explicitely by a residuum formula. We thus write 



Po,3 = Yl Po, 
ei,...,64e{±} 



3 



with an obvious definition. 

• Let us treat first pt 3 . We thus fix some C G C'nL 2 and estimate the term 

A + OVO -pi)A"(pi) 



Q f (p,p) = (2tt)- 



<P'(Pl -P2)A~(p 2 )y / (P2 -P3)A~(p 3 ) 

(^(p) + ^(P2))(E(p) + £;(p3)) 



C(P3 -P), 



by 



\Q C (p,p)\ < (2tt)- 




dpi 0£>2 dp3 



A + (p)^(p-pi)A"(pi) 



£(p + Pi) 2 /3 



E( Pi y/ 3 E( P2 )E(p) 
So if we follow the method used above, we obtain 



C(P3 ~ P) 



I Pas 



< 



3Ca/5 6 (5 6i4 ) 2 (C 6 ) 
4tt 



Wp'WI 



Now, it is easily seen that p £ ++ , p 0i3 + , pjQ" 4 " , p £ , Po,3 ++ > Po,3 
and Pq 3 ^ can be treated by exactly the same method. 
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• Let us now treat for instance p^ . Thanks to the residuum formula, we 
have to study 

Qc(p,p) = (2tt)- 6 /// dpidp 2 dpaA + (p)&(p-p 1 )A + (p 1 )y , (p 1 -p 2 )A-( P2 ) 



1 



(E(p) + E(p 2 ))(E(p) + E(p 3 ))(E( Pl ) + E( P3 )) 

If we now use the same method as above for each of the two terms of this sum, 
we arrive at 

ii ++ — ii n 3C M S6(SQ A ) 2 (C e ) 4 / 3 

llPo,3 Ik s * t- ||p Ik- 

This is easily generalized to the study of Pq 3 ++ , Po~3~ H , Pot > P03 anc ^ 

n-++- 

Po,3 

Summing now all these terms, we obtain 

ii ii / nn 3CAf5 , 6 (5' 6i 4) 2 (C' 6 ) 4 , |3 15Cm5 , 6 (S 6i 4) 2 (C 6 ) 4 .. /|l3 

llpo,3lk s ||p || £ = up || c 

47T 7T 



and 



with 



IMk < C P3 (CrV^IQIq + 2^F|p'|k) 3 , (65) 



2^ g) 3 



n „S 6 C 6 ( 5 6 C 6 V 20C M (g 6 ,4) 2 (gi 

15C A/ 5 6 (5 6 , 4 ) 2 (C 6 ) 4 



7r(4vrC ( 



,3 

coy 



□ 



4.3.5 List of constants 

We have used many constants in this proof. A summary which should help the 
reader to follow our arguments is provided in Table 1. 



Appendix: Derivation of the BDF energy 

In this section, we recall some basics about the second-quantization in no- 
photon QED and explain how the BDF energy £ is derived from this theory, 
as a mean- field approximation. We mainly follow the method of Chaix-Iracane 
[01 CU, but with the notation of [Ml El EDI- See also [HI for more details 
concerning the polarization of the vacuum. To simplify the presentation, we 
introduce P° := P° and Pi : = 1 - P°. 
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Constant 


denned in 




Lemma |3] 


Cr 


Equation (|34|) in Lemma |S] 


«i(A) 


Equation (|5*2)) in Proposition ITUl 




Equation (|56|) 


K 2j Cq 2 , C p2 


Lemma 1131 


Cm 


Lemma 1141 


Cq, Cp, C P4 


Lemma 1151 


«n, n > 4 


Lemma 1151 


Cp3' K 3 


Lemma ED 



Table 1: Constants used in the proof of Theorem 01 



Free particles, Fock space, free vacuum 

We first introduce := F%H&. and := CP^Hk which are called respec- 
tively the free electron and the free positron state subspace. C is the charge- 
conjugation operator defined by Cip = i(3a2')p- We define = J 7 ^ = C 
and J^| n) = ALiH^' ^ = AfcLi-^ for n,m > 1. The space of n free 
electrons and m free positrons is then defined by _7 r ( n > m ) = J?vy anc j th e 
associated Fock space is 

oo 

T M . (66) 

n,m=0 

For any / G TC\, the free electron (resp. positron) annihilation and creation 
operators ao(/) and ag(/) (resp. &o(/) and &o(/)) are defined as usually 
They fulfill the Canonical Anti-commutation Relations 

{ao(/),«o(5)} = {ao(f)M9)} = {&o(/), &o(<?)} = 0, (67) 

{a (f),a* (g)} = (f,P°g), {b* (f),b (g)} = (f,P°g), (68) 

where (-, •) denotes the usual scalar product of L 2 (M 3 ,C ). 

The free vacuum state Oq, a unit vector spanning jr( >°) = C, is uniquely 
characterized up to a phase factor by the properties 

||^olU = l, a (m = and 6 (/)O = 0, (69) 

for all / G Ha- 

The _/ieZc2 operator ^f(f) is defined on the Fock space J 7 by 

*(/) = ao(/) + W). 
In terms of *(/), the CAR become, for all (/,#) G H A , 

{*(/),*(<?)} = {**(/), **(g)} = 0, {*(/), **(g)} = (f,g)l, 
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Dressed particles and vacuum 

In this description, the free electrons and positrons are defined with respect to 
the projector P°, or equivalently the splitting TCa = T~t- ® 7~L\. We want now 
to change this definition and introduce the dressed electrons and positrons. To 
this end, we fix a new projector P on Ha, use again the notation P_ := P and 
P + := 1 — P, and introduce the dressed particle annihilation operators 

a P {f) : = *(P + /), b P (f) := **(P_/). (70) 

Similar formula can be given for the dressed particle creation operators a* p and 
bp. These dressed operators satisfy the same CAR as for the free operators 1)67(1 
(|68[). We also introduce the dressed electrons and positrons state subspaces 

w£ = (i - p)Ha, n p = ph a . 

Now, the main question is to know if there exists a dressed vacuum £lp in 
the Fock space T . This state has to be a solution to the analogue of 1)69)1 

||S!p||^ = 1, ap{f)n P = and b P (f)Q P = (71) 

for all / G Wa- The answer is given by the celebrated 

Theorem 4 (Shale-Stinespring Theorem |52j). There exists a dressed vac- 
uum Qp in the Fock space T satisfying f\71[ if and only if P — P° is a Hilbert- 
Schmidt operator. In this case, ftp is unique up to a phase factor. 

There are many proofs of this Theorem in the literature, see, e.g., [54 ) 136 ) 117] 
and the references in [2U] . This result explains why we assumed in the previous 
section that P - P° G S 2 (Wa)- Notice that n P can be expressed as a rotation 
of the bare vacuum in the Fock space, f2p = Uiloi U being called a Bogoliubov 
transformation. An explicit formula for can be found in a lot of papers 

[Ml ESI HZl EE EH1 [201- 

The charge of the dressed vacuum Q,p can be easily computed 2 
(Op, Qflp) = tr [P°(P - P°)P°] + tr [P°(P - P°)P-] = trpoOP ~ p °)i 



justifying our study of Section I2TT1 Here Q is the charge operator defined on T 
by [Ml Formula (10.52)] 




where (fi)i>i is an orthonormal basis of and {fi)i<—\ is an orthonormal 
basis of TL®_ . 

2 Notice that the charge of the dressed vacuum can also be easily obtained by using the 
explicit formula of flp, which immediately shows that it is an integer |54l 1371 1^01149) . 
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Second-quantized Hamiltonian 

In the physics literature, the creation and annihilation operators are defined 
differently. For instance, instead of ag(/) w hich creates an electron in the state 
Plf, the operator aj$(x) which creates an electron at x is formally used, where 
Oq{x) = Yli^i a o (fi)fi ( x )i (fi)i>i being an orthonormal basis of TL° + . The op- 
erators ap(x) and bp(x) are defined similarly. We shall now use this formalism 
which is also the one of Formally, the CAR (|67l 68|) are equivalent to 

{a (x), a (y)} = {a (x), b (y)} = W (x), b (y)} = {b (x), b (y)} = 0, (72) 



{a (x),a* (y)} = P$(x,y), {b* (x),b (y)} = P°(x,y). (73) 

We now start with writing down the formal unregularized no-photon Hamil- 
tonian 

e ur = / dx^(x) D ^ { x) +- 2 fdxf dy rm^y)m, (74 ) 

which acts on the Fock space T. As explained for instance in [Slj , the free vac- 
uum may not belong to the domain of this formally defined operator. Therefore, 
the expression (|74j) is renormalized by using a procedure which is called "normal 
ordering" , denoted by double dots : — : P o . In each product of annihilation and 
creation operators, the Oq and b$ are moved to the left as if they anticommute 
with the ao and b^. For instance 

: V*(x)9(y) : P o = a* {x)a {y) + a* (y)b {x) + b (x)a (y) - b* {y)b (x) 

= V*(x)V(y) -P°(x,y). (75) 

As a first step we thus regularize H ur as done in Chaix and Iracane [HI 
Sections 3.5 and 4.1], namely we normal order with respect to the free projector 
P°, 

at f : ty*(x)fy(x)$>*(y)fy(y) : P o 



J dx : ^*(x)D av ^(x) : P o +| / dx J dy 



i x - yi 

(76) 

This kind of regularization, which follows ideas of Dirac [131 114j . is standard 
in QED. It corresponds to the subtraction of the energy of the free Dirac sea, 
and the interaction energy with the free Dirac sea. A physical justification 
is given in [221 Section 3], on the basis of two guiding principles formulated 
by Weisskopf in |56j . The same choice is made in other studies dealing with 
vacuum polarization, for instance ^2 EH1 IS3 SHI 0H3 EH EE] • However, abetter 
choice might be possible (see the paper by Lieb and Siedentop, who propose 
another translation-invariant reference for normal ordering, in the absence of 
external field). 

Now we can express EI in terms of : — :p for some other P. Using (|75j) we 
obtain the reordering relations 

: **(z)*(y) : P o=: 9*(x)*(y) : P +Q(x, y) (77) 
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and 



: **(x)*(x)**(y)*(y) : P o = : (a:)¥*(y)tf (y) : P 

+ 2: **(x)*(x) : P T* C 4 (Q(y,y)) - 2 : **(a;)*(y) : P Q(x,y) 

+ Tr C 4 (Q(x,x))Tr C 4 (Q(y,y)) - \Q(x,y)\ 2 

where Q = P — P . Therefore we can rewrite EI with respect to an arbitrary 
dressed vacuum P formula (4.3)], 

H = / : : P * + I// ! ! ^^, 

F _ y yy F _ y| 



2 



F - 2/ 



(78) 



The last line represents the energy of the dressed vacuum P measured with 
respect to P°, whereas in the second line the vacuum polarization potentials 
appear. 



Restriction to Bogoliubov-Dirac-Fock states 

We now follow jHj and restrict ourselves to Bogoliubov-Dirac-Fock type states. 
In this approximation method, a dressed vacuum P is first chosen such that 
P — P° E &2{7i). Then a BDF state is simply a Slater determinant made 
with n electrons and m positrons defined with respect to the dressed vacuum 
P (n, m > are not fixed in this theory). This is a state of T which takes the 
form 

ip = Op(/i) • • • a*p{f n ) b* P (gi) • • • b* P (g m ) fi P , 

where (f 1} f n ) E and ...,g m ) E {H^) m are such that (/<, /_,•) = <%, 

(ft) 9j) = <5ij) and Qp is the dressed vacuum in T obtained by Theorem 0J Since 
it is easily seen that ip = Qpi where 

n m 

p' = p+i, 7=Ei/i>(/*i-Ebi><fl[7-i, 

i=l j=l 

we obtain immediately from (|78jl 

(V>|H|^) = (Op,|H|0 P /) = - P°) = £(P - P° +7). 

In formula (4.8)], this formula is expanded like in ()12l) in terms of the vacuum 
density matrix Q = P — P° and the density 7 of the dressed particles. 
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